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Summary 

The relations connecting perturbations in acoustic and entropy modes in a gas affected 
by a constant mass force, are derived. The background temperature of a gas may vary 
in the direction of an external mass force. The relations are independent on time. They 
^i \ make possible to decompose the total vector of perturbations into acoustic and non- 

■^j- ' acoustic parts uniquely at any instant. In order to do this, three quantities are required, 

according to the number of modes. In one dimension, the reference quantities may be total 
perturbations in entropy, pressure and velocity. The total energy of flow is determined. 
The examples of dismemberment of the total field into acoustic and entropy parts relate to 
the unperturbed temperature of a gas which linearly depends on the spacial co-ordinate. 
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1 Introduction 

In experiments of 1994 (Transarctic Acoustic Propagation Experiment, TAP), coherent acoustic 
transmissions via Arctic basin were studied. It was designed to monitor changes in Arctic Ocean 
temperature and sea ice properties. TAP demonstrated that the low frequency (19.6-Hz) signals 
propagated with both sufficiently low loss and high phase stability to support the coherent 
pulse compression processing and the phase detection of the signals. These yield time delay 
measurements an order of magnitude better than what is required to detect the year changes 
in travel time caused by interannual and longer term changes in Arctic Ocean temperature. 
The TAP data provided propagation loss measurements to compare with the models to be used 
for correlating modal scattering losses with sea ice properties for ice monitoring. The travel 
times measurements indicated a temperature changes in the Arctic atmosphere, which has been 



confirmed by direct measurement from icebreakers and submarines, demonstrating the utility 
of acoustic thermometry [TJ. Climatological atmospheric velocity models predict infrasound 
signals from sources that occur at mid-northern latitudes. Such infrasound data have been 
used to locate a bolid explosions in space and time. In [2] authors analyze travel time picks 
and use 3-D ray tracing to generate synthetic travel times based on various atmospheric models 
to show that the seismic network data instead reveal a predominant propagation direction. A 
sudden stratospheric warming event that reversed the zonal wind flow explains propagation 
properties. 

A time-dependent, nonlinear, fully compressible, axisymmetric, f-plane, numerical model is 
used to simulate the propagation of acoustic waves in the mesosphere and thermosphere by 
intense deep convection in the troposphere J3]- The simulations show that major convective 
storms in the tropics launch acoustic waves into the mesosphere-thermosphere directly above 
the storm centers . The principal feature of the overhead acoustic wave field in the period 
interval of 3 to 5 min is a trapped oscillation below about 80 km altitude with a period of 5 
min and a nearly vertically propagating wave with about a 3-min period above this height. 
The vertically propagating thermospheric acoustic oscillations are waves propagating upward 
from the thunderstorm source through the stratosphere- mesosphere. These predominantly 3- 
min waves are strongly driven for about 30 min after the storm event and weaken with time 
thereafter. The vertically propagating thermospheric acoustic waves may be the source of the 
F-region 3-min oscillations. Intense acoustic disturbances directly above thunderstorms may 
also be responsible for localized heating of the thermosphere [9]. 

There is an important ingredient on which a direct problem of sound propagation or inverse 
problems (such as acoustic thermometry) are based. It is a necessity the proper division of an 
atmospheric disturbance onto the characteristic types and hence a possibility to formulate a 
problem of acoustic wave generation propagation and reflection as one of them. The theoretical 
and numerical models as e.g. of [3j describing dynamics of gases and liquids affected by exter- 
nal forces are of great interest in geophysics, meteorology and wave theory [3]. The authors 
believe that the theoretical models are more desirable than the numerical methods, which are 
time-consumed and require large computer power and special converging and stable numeri- 
cal methods. Numerical methods, if necessary, are much more economical and successful, if 
they rely on the reasonable theoretical models. Moreover, some important problems may be 
resolved only theoretically. Some of these problems are considered in this study. Among other, 
they are: how many types of fluid motion exist, what they are and what energy associates 
with every type of motion. The external forces make the background of a fluid non-uniform, 
with background density, temperature and pressure depending on spatial coordinates. That 
essentially complicates the definition of linear modes (motions of infinitely small amplitude) 
taking place in such non-uniform media [5] . The number of roots of dispersion equation, if it is 
possible to determine them, or branches of possible types of motion, equals number of the con- 
servation equations [§],[I]- Even in the simplest case of flow in one dimension, the dispersion 
relations may be introduced over all wave-length range only if the background pressure and 
density depend on the coordinate exponentially [1] . Anyway, there are three types of motion 
in one dimension: two acoustic branches and, if the thermal conduction of a fluid is ignored, 
the stationary (entropy, non-wave) mode. In the flows exceeding one dimension, the buoyancy 
waves appear and interact [I], [5]. 

In this study, modes of one-dimensional flow are determined by independent on time rela- 
tions linking hydrodynamic perturbations. They are fixed for any mode and valid for arbitrary 



dependence of the background temperature on co-ordinate. These relations give possibility 
to distinguish modes analytically at any instant, to conclude about energy of everyone from 
them and to predict their dynamics. We will consider volumes of an ideal gas affected by a 
constant mass force, with variable temperature of the background. The first results allowing 
to distinguish modes due to relations of specific perturbations, were obtained relatively to the 
motion of atmospheric gas affected by gravity in [5], [8]. Fluids different from ideal gases, in- 
cluding liquids, are briefly discussed in Concluding Remarks. It is undoubtedly of importance 
in applications of meteorology and atmospheric dynamics. The case of an ideal gas with con- 
stant temperature is considered in Sec. 3. The examples relating to the linear dependence of 
the background temperature on vertical co-ordinate in the field of constant gravity force, are 
discussed in Sec. 4. 

2 Statement of problem 

2.1 Conservation equations and the total energy 

The equations governing fluid in absence of the first, second viscosity and thermal conduction 
manifest conservation of momentum, energy and mass. They are generally nonlinear. We start 
from the linearized differential equations in terms of variations of pressure and density, p' and 
p' from hydrodynamically stable stationary functions p, p, which are not longer constants: 
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The mean flow is absent, so that its velocity equals zero, V(x,y,z) = 0. The external force 
will be described by the gravity acceleration g = (0, 0, g), though it may involve readily other 
mass forces including non-inertial ones. The flow of an ideal gas is considered, which internal 
energy e in terms of pressure and density takes the form 

- " (2) 
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where 7 = C p /C v denotes the specific heats ratio. Eqs (1) describe gas motion of infinitely 
small amplitude. The background pressure and density follow from the zero order stationary 

equality, 

dp(z) 

— = ~ 9P{Z) - (3) 

The background quantities supporting the equilibrium distribution of temperature T Q (z), take 
the form 



It is well-established, that distribution of the background temperature is equilibrium if the 
parameter of static stability v(z) is positive [I], 

^)= 7 -l + 7 ^|M>0. (5) 

For some reason, it is convenient to introduce the quantity tp' instead of perturbation in density, 

p =p -7-p. (6) 
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is invariant, where 
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v = {— oo < x, y < oo, < z < /i} , (8) 

and ^i may be infinity. It readily follows from Eqs (l)-(7), that 

?1 = - /" ^ • (p^)efo = - / p'^da = 0, (9) 

dt J J a{v) 

where a is a surface circumscriptive the volume v. The invariance of e manifests conservation 
of the total energy of gas, which includes kinetic, barotropic and thermal parts. For e to be 
invariant, there is a certain freedom to establish the boundary conditions at z = and z = h: 
V z (z = 0) = V z (z = h) = 0, p' is any smooth function (impermeability condition across the 
boundaries) , or, for example, V z (z = 0) = 0, p'{z = h) = 0. 

3 Spectral problem and definition of modes 
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Introducing the new set of variables, 
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one may readily rearrange Eqs (1) into the following set: 
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where r\{z) = jMl- Eqs (11) determine the spectral problem. 



3.1 Case of the constant background temperature 

The analytical analysis of the dispersion relations and determined by them modes may be 
proceeded in the case of constant To and hence constant H . In this case, 77 = 1 and v — 7 — 1. 
The following matrix equality represents the system (11), 

d ( d d d \ 

«*<*«> = H«? W*'™ (12) 

where 

* = [U X U y U z P $] T , (13) 

f= (x,y,z) and L is the matrix operator depending on spacial partial derivatives. 

The condition of algebraic solvability of Eqs (12) may be established by the Fourier trans- 
formation using the basis functions exp(ik x x + ik y y + ik z z), \&(r, t) = J*f° exp(— iujt + ik x x + 

Det\\iul-l(k x ,k y ,k z )\\ = 0, (14) 

where I is the unit matrix, 1 represents the matrix operator L in the space of Fourier trans- 
forms. The dispersion equation (14) determines five roots of dispersion equation, or dispersion 
relations, everyone responsible for an especial type of gas motion. There are four wave modes, 
denoted by indices 1, 2, 3, 4, and the entropy mode, marked by 0: 

coo = 0, (15) 
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E = \f (p*0\ 2 + T 1 T7- + — ^. rT . ) dv, (17) 



Vectors of perturbations \l/ n (n = 0, . . .4) correspondent to eigenvalues iu n , form the Hilbert 
space L 2 (v) with the scalar product 

<*», *™> = / ( Pot n ■ f m + ^f- + , *"*" ) dv (16) 

and the invariant 
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where U , P and $ represent a sum of all parts of the eigenvectors, 
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The set \l/ n form complete set of eigenvectors. That is true for the self-adjoint boundary 
conditions. It may be readily established, that iL is symmetric in L 2 (v) (0 < n, m < 4): 

(zL#„, * m > - (* m , iLV n ) = i <b (Pjf m + P m tf n ) da = 0. (19) 

Jer(v) V ' 

The most important physically condition of impermeability at the upper and lower boundaries, 
z = and z — h, is self-adjoint: U z (z = 0) = U z (z = h) = 0. It is follow from Eqs (16)— (18) 
that E may be decomposed in to specific energies, E = J2 n =o ^n- 



4 Decomposition of acoustic and entropy modes in one- 
dimensional flow. General case 

In one dimension (k x = k y = 0) and constant To, the dispersion relations (15) determine three 
modes, or, in the other words, possible motions in the gas. Two of them are acoustic, describing 
different direction of sound propagation, and the last one is stationary, or entropy mode. It is 
responsible for stationary variations in pressure and density which result in stationary change 
in the gas temperature. In the absence of mass force, this mode would be isobaric. In the 
case of variable T , the dispersion equation, valid over all wave-length range, can not longer be 
written on, but the modes may be determined by relations linking field perturbations in this 
case as well. The completeness of the set of eigenvectors allow to represent the total vector of 
perturbations as a sum of acoustic and entropy vectors at any instant, 
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where index a denotes summary acoustic vector. Eqs (11) yield the relation connecting P(z,t) 
and $(z,t) in both acoustic branches, 
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and the link for the stationary entropy mode 
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The scalar product of eigenvectors is given by equality 
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The total vector of perturbations is a sum of orthogonal modes 
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Excluding P from the system (24), one obtains the equation describing R(z,t) = ^0r^ a (z, t): 
1 + 2tf (0)^ - 4H*(0)r?(z)£:) R{z, t) = ^ (W ( 7 - 2 - 2^(0)^)1) P 
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Its solution consists of that of the homogeneous equation and a partial solution, 

R(z,t)=R 2 f R^Diz 1 \t)dz' - i?i f R 2 (z')D(z',t)dz' + C 1 R 1 + C 2 R 2 , (26) 
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where Ri(z), R 2 (z) are determined by 
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i2(z, t) uniquely determines P a and stationary quantities $0 = & _ &a and P in accordance to 
Eqs (24). 

5 Linear dependence of H on z 

5.1 Explicit formula choice 

In this case, 77 = 1 + az, a is some non-zero constant, and v — 1 — 7 + / ~faH(0), and R\, R 2 
take the following form 



l/2aff(0) 



W Vl + a*/ U \l + a^ l + aif(0) V ; 

Eqs (25)-(29) allow to distinguish uniquely acoustic and entropy mode in any vector of total 
perturbations. Some simple conclusions follow immediately from Eqs (24). 

5.2 Contribution of only entropy mode in the total perturbation 

In this case, 

$ = % = i- 1 -^- + 7# (0)(1 + az)^\ P. (30) 

In illustrations, we take P = Pq in the form of a Gaussian impulse (a) and its derivative (b), 

(«,) P = ,exp(-(, - tf/PHQn M P = _ 2 , «P(-^-^/ o y)')^-^) , (3!) 

where /3, 7r denote the characteristic width of the impulse in units H and its amplitude. It 
would be superfluous to mention, that relations between field perturbations specifying every 
mode, Eqs (24), are valid at any instant. So, in this subsection and two subsections below, 
we do not determine the time which the samples of perturbations relate to. For definiteness, 
z = 3H(0), /3 = 0.3, and 7 = 1.4. 





a b 

Fig.l Excess pressure and entropy in the entropy mode for different aH(0) (—0.1,0,0.1' 

5.3 Contribution of only sound in the total perturbation 

In accordance to Eqs (24), 



P = Pa 
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+ 7 #(0) — I t— <I>. 
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(33) 



The values of Zq, (3, and 7 are the same as in the above subsection. 
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a b 

Fig. 2 Excess pressure and entropy in the sound for different aH(0) (—0.1,0,0.1). 

5.4 Zero total entropy 

The specific perturbations in pressure for the stationary mode is taken in both forms (a,b) as in 
Eq.(31), and the correspondent entropy is given by Eq.(30). Perturbation in entropy for sound 
is $ a = — $0; an d the excess pressure relates to $ a by means of Eq.(32). Fig. 3 shows excess 
pressure for acoustic and entropy modes in the cases (a) and (b) for different a. 
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.0 H(0) 



Fig. 3 Case of zero total entropy. aH(0) takes values (—0.1,0,0.1). 

6 Concluding remarks 

The choice of reference quantity in sound ($ a ) was made due to simplicity of relation linking 
P a and $ a for both acoustic branches. Note that in the limit of homogeneous gas (g — > 0, 
H — )■ oo, but the product gH(0) remains constant), •ygH(0) is the squared sound of velocity 
in the uniform ideal gas (about the limit look the numerical investigation [13J). In this case, 
<f' means the quantity proportional to perturbation in the entropy, ip' = (7 — l)pT s'. It is 
identically zero in both acoustic branches and is not longer suitable to be a reference quantity 
in them. Instead of, perturbation in density or in temperature may be chosen. In order to 
conclude about velocity of the acoustic mode, the knowledge of relation linking it with $ a 
is required. The relations differ in sign for different acoustic modes. That follows from the 
conservation system (11): Ui tZ (z, t) = K$i(z,t), U 2jZ (z,t) = —K§ 2 (z,t), where K is some 
integro-differential operator. It may be concluded from these equalities and Eqs (24), that if 
<J>i(z, t) = — $ 2 (z, t) and P±(z, t) = —P 2 (z, t), the total field is represented by the entropy mode 
and acoustic field with non-zero velocity U z (z, t) = Ux iZ (z, t) + U 2tZ {z, t) = Ui ;Z (z, t) = U 2>z (z, t) 
(hence, in general, the non-zero kinetic energy). Zero total kinetic energy means inevitably 
equal perturbations in pressure and entropy of both acoustic branches, P±(z,t) = P 2 (z,t), 
$>i(z,t) = $> 2 (z,t). The relations (24) are valid at any instant. In the other words, though 
perturbations in every quantity varies with time, links inside every mode remain independent 
on time. 

To make approximate evaluations of K, one may consider the waveforms which spectrum 
includes mostly large wavenumbers k z , k z H(0) » 1 [?]. That allows to expand acoustic eigen- 
values in the Taylor series and do obtain finally relations for every acoustic branch, propagating 
over the isothermal gas: 



Ui, z (z,t) 



91 



^o(7-l)(7^(0)) 3 / 2 i 
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(34) 



U 2 , z {z,t) 
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$ 2 (z',t)dz'. 



^ (7-l)(7^(0)) 3 / 2 7o 

The simple conclusion from Eqs (34) is that if j Q (&i(z',t) + $> 2 (z',t))dz', the kinetic energy at 
this instant is zero. 
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The relations linking U z with P and $ in a fluid affected by a mass force, are integro- 

differential. The exact links of excess pressure, density and velocity in unbounded volumes of 

gas with constant H, are derived exactly with regard to one-dimensional flow by one of the 

authors in [HI QUI EH] • They are valid for an ideal gas or any other fluid. For fluid different from 

an ideal gas, the equation of state, Eq.(2) should be corrected [10]. That corrects also dynamic 

equation governing an excess pressure, the second one from Eqs (1), and therefore, definition 

of modes and even definition of P,$ and U (Eqs (10)). It was shown in the study [10], that in 

the case of constant H = H(0), the relations (10) should be completed by the factor exp(£z), 

where 

AgH + B 

^~ 2H(gH-BY [ } 

where A and B follow from the leading-order series of excess internal energy, 

e! = A V = + £=. (36) 

P P 

£ equals zero only in an ideal gas. Estimations for water at normal conditions result in £ m 
0.85/2H. So that, even in a linear fluid flow over the background of constant temperature, the 
results can not be generalized by common replacing 7 by c 2 p(0)/p(0), where c denotes the sound 
velocity over the fluid without mass force, and p denotes the unperturbed internal pressure in 
it. 

Illustrations of Sec. 5 relate to aH(0) taking values —0.1, or 0.1. In the model of standard 
atmosphere, there are some extended domains of almost linear dependence of H on z: aH(0) is 
about —0.2 over the domain of z between and 10 kilometers, aH(10) ~ between 10 and 20 
kilometers and aH (30) ps 0.1 between 30 and 45 kilometers. These quantities may considerably 
vary depending on season, daytime and meteorologic conditions. 
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